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We derive a strong upper bound on the amount of Primordial Black Holes (PBHs) 
that can still be present in the Universe. Gravitational capture of PBHs by the 
Milky Way stars during their formation and subsequent accretion would produce a 
dramatic depletion of disk stars and especially of white dwarfs, unless the average 
cosmic density and mass of PBHs are severely constrained. Our finding also helps 
to discriminate among the various production mechanisms of PBHs. Moreover, we 
show that a star becomes overluminous before its disappearance into a PBH for a 
time span independent of its mass, thereby providing a characteristic observational 
signature of the considered scenario. We stress that our result allows for the existence 
of stellar-mass black holes in a mass range that is forbidden by standard stellar 
evolution. 

PACS numbers: 



The formation of black holes in the early Universe - called Primordial Black Holes (PBHs) 
- is regarded as an inescapable implication of the standard cosmological model. As recog- 
nized long ago, PBHs should be produced by the collapse of overdense regions when the 
Hubble radius equals their Schwarzschild radius [lj]. In addition, several other mechanisms 
are expected to trigger PBH formation, like the various phase-transitions occurring during 
the cosmic evolution, the collisions of cosmic loops as well as of domain walls, the infla- 
tionary reheating and a temporary softening of the cosmic equation of state ft Moreover, 
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their existence turns out to be a general feature of brane cosmology arising in theories with 
compactified extra dimensions 0]. 

A most striking feature of PBHs is the phenomenon of Hawking radiation [4j, which leads 
to their evaporation on a characteristic time scale t ev ~ 10 10 (Mpbh/10 15 g) 3 yr. Hence, PBHs 
with mass M PB h < 10 15 g should have already evaporated and the associated emission could 
have affected the cosmic evolution in several respects. Extensive work has addressed PBHs 
close to the evaporation limit Mpbh — 10 15 g, since in this case the Hawking radiation may 
produce observable effects in the present Universe [5(. So far, no positive evidence for PBHs 
has been reported and this circumstance allows to set upper bounds on their present density. 
Actually, the strongest one comes from EGRET observations of the gamma-ray background 
in the energy range 30 MeV < E < 120 GeV and implies that the contribution to the cosmic 
density parameter Q from PBHs with mass close to the evaporation limit is fipsH < 10~ 8 |6(. 
We stress that no similar constraint is available for PBHs with Mpbh > 10 15 g because they 
are still present today and their Hawking emission is negligible. 

As a matter of fact, PBHs with mass Mpbh > 10 15 g have attracted specific attention 
in the past. Since they behave as dynamically cold objects, they are very good candidates 
for the Cold Dark Matter (CDM) that dominates the mass budget of the Universe pj. For 
instance, it has been argued that PBHs in the mass range 5 ■ 10 15 g < Mpbh < 10 21 g can be 
produced by a primordial spectrum (with a characteristic mass scale) of inflationary origin 
and can make up a significant fraction of CDM [8]. A more radical option is that PBHs 
could be the dominant form of CDM in galaxies, thereby triggering their formation and so 
dispensing with the need of new elementary particles as CDM constituents 0]. Finally, the 
proposal has been put foreward that PBHs might ultimately produce the supermassive black 
holes residing in the galactic nuclei 10 . 

Surprisingly enough, little attention has been devoted so far to the astrophysical effects of 
PBHs with mass Mpbh > 10 15 g, most probably on the assumption that they hardly interact 
with the surrounding low-density matter and in particular that they undergo practically 
no accretion from the interstellar medium. Yet, if PBHs are present in sufficiently large 
number, gravitational capture by stars in the Milky Way disk during their formation process 



becomes unavoidable [111 . Il2| . Our aim is to explore the astrophysical implications of this 
phenomenon. Our conclusions set a strong constraint on the amount of PBHs in the present 
Universe, while leaving room for the existence of stellar-mass black holes in a mass range 
that is forbidden by standard stellar evolution. 

We start by remarking that much in the same way as it happens for nonbaryonic CDM 
which is currently thought to make up the galactic halos, PBHs of any mass above the 
evaporation limit are expected to undergo the same collapse process just because they are 
dynamically cold objects. We recall that in first approximation galactic dark halos are 
currently described by the Cored Spherical Isothermal (CSI) model, with density profile 
p(r) = a 2 / (a 2 + r 2 ) (a is the core radius) and Maxwellian velocity distribution with constant 
one-dimensional dispersion a [l3]. Note that a is independent of the mass of any halo 



constituent because the equilibrium distribution is achieved through violent relaxation 13 



So, it looks natural to assume that the distribution of PBHs inside the Milky Way should 
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be described by the CSI model as well, which dictates 

Pcdm(-R, z) = Pcdm(-Ro, 0) — £,2 ^ z 2 ^j i (1) 

where cylindrical Galactocentric coordinates are used and Rq ~ 8kpc denotes our Galac- 
tocentric distance. Moreover, we expect the overdensity of PBHs in the Milky Way (as 
compared to their average cosmic density) to roughly equal that of nonbaryonic CDM, 
thereby implying npBn(R, z) Mpbu — Pcbm(R,z) (^pbh/^cdm), where we are assum- 
ing that PBHs have a characteristic mass Mpbh- Various studies yield a ~ 6.4 kpc, 



Pcdm(-Ro, 0) ~ 0.3 (GeV/c 2 ) cm 3 and cr ~ 1.6 • 10 7 cms 1 as preferred values [141 ] . In 
addition, the CDM contribution to Q is ^cdm — 0.3. Hence, we get 

n PBH (/2, z) * 1.8 • 10- 39 (^1) ( f t l 2 Rl \ fip B H cm" 3 , (2) 



Mpbh; W + R 2 + z\ 
and so the distribution function of PBHs in the Milky Way is 

t PBH {R,z,v) - ( 2na2 y /2 e . (6) 

Next, let us go back to the epoch when disk stars formed in the Milky Way - which 
occurred a few Gigayears ago - and let us focus our attention on a single protostar cloud 
located at (R, 0), i.e. a Giant Molecular Cloud (GMC) of typical mass M GMC ~ 10 39 g 
and linear size -Rgmc — 10 pc. According to current wisdom, it produces many stars upon 
fragmentation associated with the isothermal gravitational collapse. Therefore, the number 
of unbound PBHs initially contained in such a GMC is simply its volume times the number 
density given by Eq. (jSJ), namely 

W ,o). 2 ,.io-(i£5)(J±||) fipBH . (4, 

Owing to the GMC collapse, a PBH inside it feels a time- decreasing gravitational potential. 
So, its en ergy gets lowered, thereby allowing a fraction of PBHs in the collapsing GMC to get 



captured [15|. Below, we proceed to a quantitative estimate of this effect, following closely 
the analysis of Steigman, Sarazin, Quintana and Faulkner (SSQF) 0. As stressed by these 
authors, the realistic situation is the one in which the collapse time t c = $gmc/(^gmc/^) 
is much longer than the free-fall time tg = (R^ MC /GM GMC ) l l 2 , where $qmc denotes the 
cloud gravitational potential. Assuming that the PBH distribution function in the considered 
GMC is given by Eq. (EL by going through the same steps of SSQF we find that the total 



number of PBHs captured in the course of the GMC collapse is 17 



«(^0) * 1.1 i^^) 3/2 N^,(R,0) , (5) 



that is 
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Since even the isothermal collapse has an efficiency of at most 40% 18j, we expect that 
no more than 10 5 stars should form out of a single GMC As a consequence, the average 
number of PBH gravitationally bound to each such star should be 

^(«,0)-1.9.10^(i£i)(J±|) npB „, (7) 

where the dilution factor 5 accounts for the possibility that some PBHs end up in the 
intracluster gas rather than in a star. Thus, we see that on average every star in the Milky 
Way disk contains at least one PBH provided that 

Mpbh < 1-9 ■ 10 27 5 + fipBH § ' ( 8 ) 

where now R stands for the star's Galactocentrc distance. 

Once a PBH with mass Mpbh gets captured by a star, it sinks towards the center and 
starts accreting. In such a situation, its mass increases with time and will be denoted by 
Mpbh(^)- Given that PBHs are produced in the early Universe and since we are dealing with 
a stellar population formed a few Gigayears ago, we can assume Mpbh — Mpbh(O), while the 
present value of a captured PBH mass is Mpbh(^h), where tn ~ 1.3 • 10 10 yr is the Hubble 
time. Similarly, the resulting accretion luminosity L acc (t) increases, eventually becoming 
comparable to the original stellar luminosity L* - say, L acc ~ 0.1L* - on a characteristic 
time scale t\ um . It may also happen that the whole star gets swallowed by the PBH on a 
characteristic time t sw . Below, this picture will be worked out in a quantitative fashion. 

As a preliminary step, we recall a few basic facts that are instrumental to our subsequent 
analysis. The PBH Schwarzschild radius is R s ~ 1.5 ■ 10~ 13 (Mpbh/10 15 g) cm, its Hawking 
temperature is Th — 1.2 • 10 11 (10 15 g/MppH) K and the corresponding luminosity is Lh — 
3.3- 10 15 (10 15 g/MpBH) 2 ergs -1 . Although Th can largely exceed that of stellar matter, Lh is 
totally negligible in comparison to L,. Within the Bond, accretion theory 0, the accretion 
radius of any object of mass M is i? acc ~ 1.5 • 10~ 4 (M/10 15 g) (10 4 K/T b d) cm, while the 
accretion rate is given by 

where pbd and T bd denote the density and the temperature, respectively, of matter at -R acc . 
Eq. (j5J) can be trivially solved to yield 

The associated accretion luminosity is 

2 dM , , 

^acc ^ CC — , (11) 



where e denotes the efficiency of the accretion process. The issue of the efficiency in spherical 
accretion has been widely debated in the past and e was found to depend on the accretion 
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rate [20J. In the regime of large optical depth - typical of the stellar interior - it turns out 
that KT 5 < e < 1(T 4 . Whenever necessary, we will conservatively take e ~ 10 5 . Generally 
speaking, as the accretion proceeds L acc increases until the Eddington luminosity 

Le ~ 6 ' 5 ' 1019 (l^i) 6rgS_1 (12) 

is attained, after which the accretion process gets Eddington-limited and the accretion rate 
becomes 

dM „ „ , „ R _! / M \ _ x 



a - 2 - 310 e U^j gyr ■ (13) 

It proves convenient to represent the solution of Eq. ffTB"]) in the form 

M(t) = M(0) 10 x(t) , (14) 

with X(t) ~ 1 ■ lO^e- 1 (t/yr). 

Coming back to our main line of development, we inquire about the fate of a disk star 
that has indeed eaptuted a PBH 0. Specifically, we will address both the case of a Sun-like 
star and that of such a star having evolved into a white dwarf. 

As far as the behaviour of a Sun-like star is concerned, we assume for defmiteness 
M* ~ 2 • 10 33 g, L* ~ 4 • 10 33 ergs~\ p bd ~ lOgcm" 3 and T bd ~ 10 7 K. Then 
the PBH settles at the centre because of dynamical friction on a time scale t% ~ 
7.4 • 10 9 (fpBH/10 7 cms -1 ) 3 (10 15 g/M PBH ) yr [131 ]. the accretion radius is _R acc ~ 1.5 ■ 
1CT 7 (MpBH/10 15 g) cm [22j and the accretion luminosity turns out to be 

Li^8.6.10» f (^f|) 2 erg 8 -'. (15) 

The comparison of Lf cc (t) with L^{t) shows that three possibilities can be realized. (1) 
Accretion is never Eddington-limited: this happens for L^ c (t H ) < -^e(^h), namely for 
MpbhM < 7.6 • 10 19 e~ l g, i.e. M P bh(*h) < 7.6 ■ 10 24 g. Owing to Eq. ([TO]), such a 
condition translates into M PB h < 2.6 ■ 10 18 g (regardless of e). Correspondingly, the presence 
of the PBH inside the star is totally harmless, since L^ cc (tu) <^ L* and Mpbh(^h) ^ M*. 
(2) Accretion is always Eddington-limited: this takes place for L s ^ cc (0) > Le(0), namely 
for Mpbh > 7.6 • lO^e-ig, i.e. M PBH > 7.6 • 10 24 g. Clearly, t sw is fixed by the con- 
dition M PBR (t sw ) ~ M*. Hence, by setting M(0) = M PBH = 10 y g and M{t) = Af», 
Eq. (P3J) entails t sw ~ (33.3 - F) ■ 10 9 eyr, i.e. t sw ~ (33.3 - Y) ■ 10 4 yr. Since now 
Y > 19.9 - loge, we deduce t sw < (13.4 + loge) • 10 9 eyr, i.e. £ sw < 8.4 • 10 4 yr. More- 
over, t\ nm is fixed by the requirement Lg^ium) — 0.1 L*. Accordingly, Eq. f|T2l yields 
MpBH(*ium) - 6.2 • 10 27 g and using again Eq. flU]) we find t lum ~ (27.8 - Y) ■ 10 9 eyr, i.e. 
t\ nm ~ (27.8 — Y) ■ 10 4 yr. So, for a time span of about 5.5 ■ 10 9 eyr, i.e. 5.5 ■ 10 4 yr the star 
becomes overluminous before disappearing (regardless of M*). (3) Accretion starts in the 
Bondi regime and next gets Eddington-limited for 2.6 • 10 18 g < Mpbh < 7.6 ■ 10 19 e _1 g, 
i.e. 2.6 • 10 18 g < Mpbh < 7.6 ■ 10 24 g. The turnover occurs at time t to such that 
MpBH(tto) - 7.6 • 10 19 e~ 1 g and can be found from Eq. ([Hty by setting M(0) = M PBH 
and M(t) = M PBH (t to ). We obtain t to ~ 3.3 • 10 13 [(10 15 g/M PBH ) - 1-3 ■ 10~ 5 e] yr and hence 
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we have t to < 1.3 • 10 10 yr. Accretion next gets Eddington-limited and proceeds much more 
rapidly, as described above. 

We next carry out a similar investigation for a white dwarf, taking for definiteness M* ~ 
2 • 10 33 g, L* ~ 7 • 10 30 ergs _1 , p bd ~ 10 6 g cm" 3 and T bd ~ 10 7 K. Correspondingly, i? acc 



is still given by the previous expression 22], whereas now the PBH settles at the centre 
because of dynamical friction on a time scale tjf ~ 10 -5 tf f and the accretion luminosity 
becomes L" c ^ ~ 10 5 L s ^ cc . Proceeding as above - namely comparing L%£,(t) with L-&{t) - the 
following scenario emerges. (1) At variance with the previous case, a pure Bondi accretion 
regime is ruled out, because now pu is larger by a factor 10 5 and so it would demand 
Mpbh < 2.6 • 10 13 g, which is against our main assumption. (2) The same argument implies 
that accretion is always Eddington-limited for Mpbh > 7.6-10 14 e -1 g, i.e. Mpbh > 7.6-10 19 g. 
Since M* is unchanged, t sw is the same as before. However, we presently have Y > 14.9— log e, 
which entails t sw < (18.4 + log e) ■ 10 9 eyr, i.e. t sw < 1.4-10 5 yr. Finally, we evaluate ti um from 
the condition L E (ti um ) — 0.1 L*. Correspondingly, Eq. (112]) gives MpBH^ium) — 11 • 10 25 g 
and from Eq. (1T4"1) we get t\ um ~ (25 — Y) ■ 10 9 eyr, i.e. t\ um ~ (25 — Y) • 10 4 yr. Hence, also 
a white dwarf becomes overluminous for a time span of about 8.3 • 10 9 eyr, i.e. 8.3 • 10 4 yr 
(regardless of M*) before disappearing. (3) Accretion starts in the Bondi regime and next 
gets Eddington-limited for 10 15 g < M PBH < 7.6 ■ 10 14 e" 1 g, i.e. 10 15 g < M PBH < 7.6 • 10 19 g. 
The turnover takes place at time t to such that MpBH^to) — 7.6 ■ 10 14 e _1 g and can be 
computed from Eq. by setting M(0) = M PBH and M(t) = M PBU (t to ). We find t to ~ 
3.3 • 10 8 [(10 15 g/MpBH) — 1.3 e] yr and consequently t to < 3.3 • 10 8 yr. Thereafter, accretion 
becomes Eddington-limited and proceeds much faster, as explained above. 

Our findings set a strong upper bound on the amount of PBHs in the present Universe. 
To see how this comes about, suppose that every star in the Milky Way disk captures a 
PBH during its formation process 23] . We have shown that in such a situation the star gets 
swallowed by a PBH with Mpbh > 10 20 g on a time scale much shorter than tn 2J, [25j. This 



would not be the case for lighter PBHs, since they have not enough time to engorge the star. 
Nevertheless, a similar situation would inevitably occur on a time scale much shorter than 
- for any value of Mpbh above the evaporation limit - once the original star has evolved into 
a white dwarf, owing to the enhanced accretion triggered by the increased density. In order 
to avoid such a catastrophic conclusion, we have to demand that only a small fraction of 
disk stars in the Milky Way do capture a PBH. This requirement implies that condition (jEJ) 
should be grossly violated. As we said, no more than 40 % of the mass of a GMC goes into 
stars and it looks reasonable to expect a star-formation efficiency of order 10 %, in which 
case the dilution factor is 5 ~ 0.1. Thus, we conclude that only PBHs with M PB h ^ 10 25 g 
can give a substantial contribution to the dark matter in the Universe. More generally, this 
bound strongly constraints the initial mass function of PBHs and so it helps to discriminate 
among their various production mechanisms. 

Our results also give rise to an exciting possibility concerning stellar-mass black holes. 
According to the standard theory of stellar evolution, they are the remnant of supernova 
explosions of stars more massive than about 20 M 26] and their mass invariably turns out 
to exceed 3 M . Therefore, stellar-mass black holes with mass M < 1 M Q can be neither the 
final stage of stellar evolution nor the product of neutron star-black hole transition through 
accretion. However, a PBH that swallows a white dwarf ends up with a mass equal to that 
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of the white dwarf itself, namely in the range 0.2 M & — 1.4 M & . Thus, observational evidence 
for black holes with mass M < 1 M & would point to the existence of PBHs more massive 
than 10 15 g. 

Obviously, several aspects of the scenario outlined above require further investigation. 
For instance, it is not clear to us whether an explosive event is produced when the star 
gets swallowed by a PBH. Still, the characteristic overluminous phase that precedes the 
disappearance of the star - whose time span is independent of the stellar mass - should be 
actually detectable. 

We thank an anonymous referee for criticism which has led to a substantial improvement 
of this Letter. One of us (M. R.) would like to thank the Dipartimento di Fisica Nucleare e 
Teorica, Universita di Pavia, for support. 



[1] Ya. B. Zel'dovich and I. D. Novikov, Sov. Astron. 10, 602 (1967); S. Hawking, Mon. Not. R. 
Astron. Soc. 152, 75 (1971); B. J. Carr and S. W. Hawking, Mon. Not. R. Astron. Soc. 168, 
399 (1974). 

[2] For a review, see: B. Carr, in Inflating Horizon of Particle Astrophysics and Cosmology 

(Universal Academy Press, 2005). 
[3] R. Guedens, D. Clancy and A. R. Liddle, Phys. Rev. D 66, 043513 (2002); 66, 083509 (2002). 
[4] S. W. Hawking, Nature 248, 30 (1974); Comm. Mat. Phys. 43, 199 (1975). 
[5] For a review, see: F. Halzen, E. Zas, J. H. MacGibbon and T. C. Weekes, Nature 353, 807 

(1991). For an updated list of references, see: M. Schroedter et al, arXiv:0812.0546 (2008). 
[6] D. N. Page and S. W. Hawking, Astrophys. J. 206, 1 (1976); B. J. Carr and J. H. MacGibbon, 

Phys. Rep. 307, 141 (1998). 
[7] For a review, see: B. J. Carr, Ann. Rev. Astron. Astrophys. 32, 531 (1994). 
[8] D. Blais, T. Bringmann, C. Kiefer and D. Polarski, Phys. Rev. D 67, 024024 (2003). 
[9] N. Afshordi, P. McDonald and D. N. Spergel, Astrophys. J. 594, L71 (2003). 
[10] N. Dutchting, Phys. Rev. D 70, 064015 (2004). M. Yu. Khlopov, S. G. Rubin and A. S. 

Sakharov, Astropart. Phys. 23, 265 (2005). 
[11] C. Bambi et al. arXiv:0812. 05851 (2008). 
[12] M. A. Abramowicz et al. arXiv:0810.340v2 (2009). 

[13] J. Binney and S. Tremaine, Galactic Dynamics (Princeton University Press, Princeton, 2008). 
[14] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rep. 267, 195 (1996); M. Kamionkowski 

and A. Kinkhabwala, Phys. Rev. D 57, 3256 (1998); A. M. Green, Phys. Rev. D 63, 043005 

(2001). 

[15] In order to avoid any confusion, we stress the captured is understood throughout as a synonym 

of gravitationally bound, since evaporation is totally irrelevant [l6| . 
[16] G. Steigman, C. L. Sarazin, H. Quintana and J. Faulkner, Astron. J. 83, 1050 (1978). 
[17] We note in passing that the numerical coefficient in Eq. (13) of SSQF should read 1.1 rather 

than 0.51. 

[18] J. M. Scalo, in Protostars and Planets II, ed. by D. C. Black and M. S. Mathews (University 

of Arizona Press, Tucson, 1985). 
[19] S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars (Wiley, 

New York, 1983). A. Treves, L. Maraschi and M. Abramowicz, Pub. Ast. Soc. Pac. 100, 427 



8 



(1988). 

[20] L. Nobili, R. Turolla and L. Zampieri, Astrophys. J. 383, 250 (1991). 

[21] This issue has been discussed in a different context in: D. D. Clayton, M. J. Newman and 
R. J. Talbot, Astrophys. J. 201, 489 (1975); D. Markovic, Mon. Not. R. Astron. Soc. 277, 
25 (1995); M.C. Begelman, E. M. Rossi and P. J. Armitage, Mon. Not. R. Astron. Soc. 387, 
1649 (2008). 

[22] The accretion radius is so small that taking the central values for the density and the tem- 
perature is absolutely correct for a PBH in the central region of the star. 
[23] Subsequent PBH capture by a star turns out to be negligible, as pointed out in ref. 



12 . 



[24] This follows from our previous estimate of ito and shows that the PBH capture and its settling 
at the star centre indeed occur early enough on the cosmic time scale. 

[25] This conclusion was overlooked in D. D. Clayton, M. J. Newman and R. J. Talbot, Astrophys. 
J. 201, 489 (1975). 

[26] As usual M & ~ 2 • 10 33 g denotes the solar mass. 



